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INTRODUCTION 


The  following  annual  progress  report  documents  work  at  the  University 
of  Minnesota  Heat  Transfer  Laboratory,  under  AFOSR  sponsorship,  on  the  topic 
of  heat  transfer  from  gas  turbine  airfoil  and  end-wall  surfaces.  The 
research  was  divided  into  the  following  three  subtopics:  Curvature  Effects 
and  End-Wall  and  Near-End-Wall  Heat  Transfer.  Progress  on  each  is  discussed 
separately. 

A.  The  Effect  of  Convex  Curvature  on  Heat  Transfer  and  Hydrodynamics  of 
Turbulent  Boundary  Layers 

Three  experiments  were  conducted  to  investigate  the  effect  of 
streamwise  curvature  over  the  range  from  6/R  -  0.01  to  0.0*1  and  to  study  the 
effect  of  free-stream  turbulence  on  fully  turbulent  curved  boundary  layers. 
An  attached  fully  turbulent  boundary  layer,  grown  on  a  flat  plate  is 
introduced  to  a  convex  wall  of  constant  radius  of  curvature  followed  by  a 
flat  recovery  wall  (Figure  1).  Descriptors  of  the  three  experiments  are 
listed  in  Table  1 . 

1 .  Objectives 

The  effects  of  varying  the  radius  of  curvature  on  curved, 
two-dimensional  boundary  layers,  including  the  recovery  process,  were 
investigated  to  further  understand  the  transport  mechanism  and  to  extend  the 
data  base  for  prediction  model  development.  Two  different  free-stream 
turbulence  intensity  cases  were  employed  to  study  the  effects  of  free-stream 
turbulence  intensity  on  the  curved  boundary  layer.  The  higher  turbulence 
intensity  case,  Case  3.  showed  indications  of  interactions  between  the  outer 
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(unheated)  and  the  inner  (heated)  shear  layers  for  x  >  150  cm  in  the 
recovery  station.  Significant  additions  to  the  previous  data  base  are: 

(1)  The  recovery  process,  with  heat  transfer,  was  investigated  for 
cases  of  milder  curvature  than  that  of  the  previous  studies. 

(2)  The  radius  of  curvature  was  varied  from  case  to  case  in  the  same 
study. 

(3)  Curvature  effects  of  6/R  between  those  of  mild  and  strong 
curvature  were  investigated  in  the  same  facility. 

(M)  Free-stream  turbulence  intensity  effects  on  the  curved  boundary 
layer  were  investigated  for  the  first  time. 


2.  Progress  and  Accomplishments 

The  important  observations  from  the  three  cases  of  experiments  are 
summarized  based  on  You  (1986);  You,  Simon  and  Kim  (1986a)  and  You,  Simon 
and  Kim  (1986b)  —  the  second  and  third  references  are  attached: 

(A)  General  characteristics  of  convex-curved  boundary  layers:  CASE  1 

(1)  Profiles  of  mean  velocity  and  temperature  show  shortened 
log-linear  regions  and  enhanced  wake  regions  in  the  curve. 

The  recovery  of  these  mean  profiles  in  the  log-linear  region 
is  complete  by  about  80  cm  of  recovery  length.  The  wake  is 
decreasing  slightly  near  the  end  of  the  recovery,  indicating 
incomplete  recovery. 

(2)  The  curvature  effect  increases  the  turbulent  Prandtl  number 
in  the  log-linear  region  by  about  25%  for  Case  1. 

(3)  Reduced  rates  of  growth  of  momentum  thickness  and  increased 
shape  factor  within  the  curve  due  to  the  effect  of  a  decrease 
in  the  turbulent  transport  are  observed.  The  recovery  of  the 
shape  factor  H  is  not  complete  at  the  end  of  the  recovery 
section. 

(M)  Rapid  decreases  of  St  and  Cf/2  at  the  beginning  of  the  curve 
are  observed,  followed  by  slow  decreases  within  the  curve 
(Figure  2).  Decreases  from  expected  flat-wall  values  are  20% 
for  Case  1  (6/R  «  0.03).  The  recovery  process  is  extremely 
slow  with  St  recovery  leading  Cp/2--an  increase  of  2St/Cf. 

(5)  Streamwise-normal  turbulence  intensity  profiles  (Figure  3) 
show  a  fast  response  to  the  introduction  of  curvature  with  a 
drastic  reduction  of  the  near-wall  peak.  The  self-similar 
shape  is  attained  by  station  and  continues  through  the 
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■  remainder  of  the  curve.  The  recovery  process  begins  with  an 

(increase  of  turbulence  intensity  near  the  wall  at  y/6  = 

0.1 -0.2  propagating  slowly  to  the  outer  layers.  Shear  stress 
profiles  respond  quickly  to  curvature,  reducing  to  about  55? 
of  the  upstream  flat-plate  values  for  6/R  =  0.03.  Recovery 
!  occurs  slowly,  with  the  wake  recovering  faster  than  the 

near-wall  flow. 


(B)  Effect  of  radius  of  curvature:  Comparison  between  CASE  1  and 
CASE  2 
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(1)  Though  both  Cases  1  and  2  show  the  same  trends,  the  effects 
of  the  different  strengths  of  curvature  are  clearly 
observable. 

(2)  Decreased  St  and  Cf/2  of  about  10?  from  expected  flat-wall 
values  are  observed  for  Case  2  (6/R  =  0.013)  whereas 
decreases  of  about  20?  are  observed  for  Case  1  (6/R  =  0.03) 

— Figure  2.  The  ratio  of  boundary  layer  thickness  to  the 
radius  of  curvature  is  the  appropriate  scaling  parameter  over 
the  range  from  6/R  =  0.01  to  6/R  *  0 . 04 . 

(3)  The  turbulent  Prandtl  number  in  the  log-linear  region  shows  a 
20?  increase  for  Case  2. 

(*0  "Asymptotic"  behavior,  suggested  by  Gillis  and  Johnston 
(1983)  is  observed. 

(5)  Shear  stress  profiles,  plotted  in  -u'v'/UpW2  vs.  y/6,  show 
the  realization  of  a  self-similar  shape  by  station  3  even  in 
the  weaker  curvature  case. 


(C)  Effect  of  free-stream  turbulence  intensity:  Comparison  between 

CASE  1  and  CASE  3 

(1)  The  curvature  effect  dominates  the  free-stream  turbulence 
effect  for  the  cases  studied  (6/R  =  0.03). 

(2)  In  Figure  4 ,  flat-plate  data  at  x  •  -29.9  cm  show  reduced 
wake  strength  and  about  5?  increases  in  St  and  Cf/2  due  to 
the  higher  free-stream  turbulence  intensity  (1.85?). 

(3)  Higher  free-stream  turbulence  intensity  causes  a  stronger 
curvature  effect  on  hydrodynamic  transport  than  on  thermal 
transport  (Figure  *0:  Cf/2  within  the  convex-curved  section 
decreases  more  than  St  decreases,  which  is  consistent  with 
the  decrease  of  log-linear  region  Pr^  near  the  end  of  the 
curved  section. 
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curvature,  asymptotic  turbulent  shear  stress  profiles  and 
slow  recovery. 

(5)  Similar  streamwise  turbulence  intensity  profiles  for  Cases  1 
and  3  near  the  end  of  the  curve  were  found  in  spite  of  the 
different  streamwise  turbulence  intensity  profiles  at  station 
1  (the  difference  in  profiles  was  a  result  of  the  different 
free-stream  turbulence  intensities). 

(6)  Reynolds  shear  stress  profiles  are  dramatically  influenced  by 
curvature--the  Reynolds  shear  stress  data  of  y/6  >  0.7 
becomes  negative  at  station  if.  A  self-similar  shear  stress 
profile  in  -u'v'/UpW2  vs.  y/6  coordinates  in  the  curved 
section  was  realized  inside  50?  of  the  shear  layer  thickness. 


( D)  The  curvature  effect  on  the  wake 

(1)  Coles  wake  law  is  valid  in  the  convex-curved  boundary  layer 
when  an  appropriate  wake  parameter  is  found,  using: 


II  =  U/2)  -±-  -  { ( 1  /< )  ln(  6U  /v)  ♦  C}] 

Ut  1+5/R  t 


a 


(2)  The  variation  of  the  wake  parameter  in  the  curved  flow  is 
different  than  that  observed  for  flat-plate  equilibrium  cases 
(Figure  5);  but,  the  variation  is  consistent  for  the  two 
low-turbulence  cases. 

(3)  The  effect  of  curvature  on  the  wake  parameter  shows 
consistent  trends  with  the  other  measurements:  a  stronger 
effect  for  stronger  curvature  and  an  extremely  slow  recovery. 

(*J)  The  distribution  of  the  wake  components  in  the  higher 

turbulence  case  (Case  3)  did  not  indicate  similarity  with 
Coles  wake  law. 


3.  Expected  Activities 


The  convex  boundary  layer  study  has  been  completed. 

B.  Modeling  of  Curved  Turbulent  and  Transitional  Flows 
1 .  Objectives 

(1)  Test  transition  models  against  recent  data. 
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(2)  Test  curvature  corrections  to  turbulent  models  against  recent 
data . 

(3)  Improve  transition  and  turbulent  models. 

2.  Progress  and  Accomplishments 

( 1 )  Test  of  an  all-laminar  boundary  layer  flow. 

An  all-laminar  flow  was  simulated  using  STAN5  (Crawford  and  Kays, 

1976),  a  computer  code  for  solving  two-dimensional  boundary  layer  flows. 

The  data  was  of  the  qualification  test  for  the  transitional  boundary  layer 
experiment  of  Wang,  Simon  and  Buddhavarapu  (1985).  Measurements  were  taker, 
in  an  all-laminar  flow  over  a  flat  plate  with  a  constant  heat  flux.  The 
predictions  showed  the  effect  streamwise  acceleration  on  the  skin  friction 
coefficient  and  Stanton  number. 

(2)  Modifications  to  STANS. 

STANS  was  modified  mainly  to  employ  the  transition  models  and  curvature 
modifications  to  the  ML'ri  turbulence  model.  In  the  m3in  program,  four 
transition  start  models  and  three  transition  length  models  were  added. 

Three  transition  path  models  and  two  curvature  models  were  employed  in  one 
subroutine.  The  input  programming  was  modified  so  that  the  models  can  be 
selected  and  so  that  the  turbulent  Prandtl  number  could  be  input  as  a 
function  of  streamwise  distance. 

In  addition,  the  property  subroutine  was  revised  so  that  SI  units  can 
be  used,  and  the  calculation  and  printout  of  Reynolds  shear  stress  and 
turbulent  heat  flux  was  added. 

( 3)  Test  of  transition  model. 

(i)  Transition  start  model  (Figure  6  a,  b) 

Transition  models  were  tested  against  recent  data  taken  by  Wang,  Simon 


and  Buddhavarapu  (1985)  for  freestream  turbulence  intensity  of  0.68?  and 


2.0$  over  a  flat  plate  with  a  constant  heat  flux.  In  both  cases,  the 


Abu-Ghannair.  and  Shaw  (1980)  model  predicts  the  transition  start  very  well. 
Seyb  (1972),  Dunham  (1972)  and  van  Driest  sr.d  Blumer  (1  963)  models  predicted 
a  late  start.  Abu-Ghannam  and  Shaw  was  therefore  selected  as  the  transition 
start  model. 

(ii)  Transition  length  model  (Figure  7  a,  b) 

With  the  Abu-Ghannam  and  Shaw  ( 1  930 )  start  model  and  the  Dhawan  and 
Narasimha  (1958)  path  model,  three  transition  length  models  were  tested. 

The  results  showed  that  the  Chen  and  Thyson  (1971)  and  Debruge  (1970)  models 
predict  short  transition  lengths  resulting  in  a  steeper  increase  in  Stanton 
number  in  the  transition  region.  The  Dhawan  and  Narasimha  model  yields  a 
steeper  response  than  the  experimental  data  would  indicate  but  is  superior 
to  the  other  two. 

(iii)  Transition  path  model  (Figure  8  a,  b) 

Among  the  three  path  ( intermittency )  models;  Dhawan  and  Narasimha 
(1958),  Abu-Ghannam  and  Shaw  (1980),  and  Chen  and  Thyson  (1971),  the  Dhawan 
and  Narasimha  model  gives  a  better  comparison  to  the  data  than  do  the 
others . 

(*0  Test  of  curvature  model.  (Figure  9  a,  b,  c) 

The  Adams  and  Johnston  (1983)  and  the  Bradshaw  (1973)  curvature 
modifications  to  the  MLH  turbulence  model  were  tested  against  recent 
turbulent  data  over  a  convex-curved  wall  measured  by  You  (1986)  and  You, 
Simon  and  Kim  (1986a,  1986b).  In  the  three  cases  tested,  the  Adams  and 
Johnston  model  shows  good  agreement  with  the  experimental  data. 

(5)  Test  of  the  combined  model  against  curved-wall  transition  data 
The  combined  model  with  the  three  transition  models,  Abu-Ghannam  and 


Shaw  (start),  Dhawan  and  Narasimha  (length),  Dhawan  and  Narasimha  (path), 
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and  Adams  and  Johnston  (turbulent  curvature  model)  was  tested  against  the 
convex-curved  wall  transition  data  (Wang  and  Simon,  1985).  In  the  higher 
turbulence  intensity  (2$)  case,  the  combined  models'  results  compared  well 
as  can  be  seen  in  Figures  10  and  11. 

Through  the  tests  of  the  transition  and  curvature  models  using  the 
mixing  length  hypothesis  turbulence  model,  it  was  found  that: 

(i)  the  start  of  transition  is  predicted  well  by  the  Abu-Ghannam  and 
Shaw  model  (1980) . 

(ii)  the  Dhawan  and  Narasimha  (1958)  models  predict  the  transition 
length  and  path  better  than  the  other  models  tested. 

(iii)  the  predicted  onset  of  transition  is  late  and  the  transition  path 
is  rather  rapid  compared  to  the  data. 

(iv)  the  combined  model  of  the  transition  and  turbulence  curvature 
models  is  applicable  and  accurately  calculates  the  curved  wall 
transition  data. 

3.  Expected  Activities 

The  following  will  be  continued: 

(1)  Test  of  the  selected  models  against  other  data  sets 

(2)  Improvement  of  the  models 

(3)  A  study  of  the  low  Reynolds  number  version  K-e  model 

C .  Measurement  of  the  Turbulent  Transport  of  Heat 

1 .  Objectives 

The  objective  for  this  portion  of  the  work  was  to  develop  a  probe  which 
could  be  used  to  measure  the  turbulent  transport  of  heat,  v't'. 

2.  Progress  and  Accomplishments 

The  probe  design,  fabrication  and  qualification  is  now  complete.  The 
following  summarizes. 
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Design  of  Probe.  In  designing  special  purpose  hot-wire  arrays , 
consideration  must  be  given  to  a  multitude  of  factors,  some  of  which  are 
spatial  averaging,  prong  and  shaft  interference,  end  conduction,  sensor 
cross-talk,  and  survivability.  Ideally,  one  would  like  a  probe  whose  active 
length  was  smaller  than  the  smallest  scales  of  interest  while  maintaining 
practicality  of  assembly  and  sensor  strength.  Since  the  Kolmogorov  scale  of 
laboratory  turbulence  typically  ranges  from  5  pm  to  100  pm,  sensors  of  this 
type  are,  as  of  yet,  impossible  to  build.  Compromises  must  therefore  be 
made . 

The  choice  of  wire  length  and  diameter  was  guided  by  the  four  criteria 
presented  by  Blair  and  Bennett  (19814):  (1)  Reduction  of  mean  flow  gradients. 

The  gradients  of  velocity  and  temperature  along  the  wire  were  reduced  by 
choosing  the  smallest  wire  length  possible  subject  to  other  considerations. 
The  ratio  between  boundary  layer  thickness  and  wire  height  was  a  minimum  of 
60  at  all  times.  (2)  Reduction  of  end  effects.  End  effects  were  reduced  by 
choosing  a  length  to  diameter  (1/d)  of  200,  as  recommended  by  Champagne 
(1967).  This  also  enabled  the  correction  equations  recommended  by  Champagne 
(1967)  to  be  employed.  The  active  length  consisted  of  only  the  central 
one-third  of  the  wire,  the  ends  of  the  wire  being  plated  to  reduce  end 
losses  further.  (3)  Practicality  of  fabrication  and  sensor  survivability. 

The  smallest  wires  commonly  used  for  hot  wire  sensors  are  of  2.5  pm  plated 
tungsten.  These  wires  combine  the  strength  of  tungsten  with  the 
non-reactive  nature  of  platinum.  The  wires  can  therefore  be  operated  at 
high  temperatures  with  good  survi vaoility .  Furthermore,  platinum  presents  a 
smooth  surface  to  the  flow,  minimizing  calibration  drift  due  to  accumulation 
of  dust  on  the  wire.  (M)  Assurance  of  adequate  spacial  resolution.  The 
spatial  resolution  of  the  wire  may  be  calculated  from  the  curves  presented 


by  Wyngaard  0968).  It  is  found  that  for  turbulence  scales  greater  than  10 
Kolmogorov  scales,  signal  attenuation  is  negligible. 

The  distance  between  the  prongs  is  determined  by  the  practicality  of 
fabrication,  interference  effects,  and  spatial  resolution.  The  transverse 
correlation  coefficient  falls  rapidly  with  sensor  separation,  as  pointed  out 
by  Sandborn  (1972)  and  Fulachier  (1978).  The  sensors  should  therefore  be 
placed  as  close  as  possible,  but  not  so  close  as  to  interfere  with  one 
another.  A  transverse  spacing  between  the  X-array  of  0.78  mm  was  chosen. 

The  configuration  of  the  wires  is  as  shown  on  Figure  12.  The  probe 
consists  of  three  wires  located  in  three  parallel  planes.  The  two  outer 
wires  are  orthogonal,  forming  an  X-array,  while  the  third  center  wire  is 
parallel  to  one  of  the  outer  wires.  The  two  outer  wires  are  operated  at  the 
same  overheat  (0H=1.9),  while  the  center  wire  is  operated  at  a  much  lower 
overheat  (0H=1.035).  The  instantaneous  temperature  is  then  found  from  the 
two  parallel  wires,  and  the  instantaneous  velocity  from  the  two  outer  wires, 
as  will  be  shown  in  the  data  reduction  section. 
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incorporated  into  the  response  equation  using  a  least-squares  fit  to  the 


data.  The  final  response  equation  is  given  by 


U  0  •  - 3  5  =  ;  rr.  <  7  +  h » }  T 0  •  76 


.  +  bB)w2 

(Tw  -  7-r)  Var 


For  isothermal  flows,  i.e.,  Inconstant,  the  response  equation  reduces  to 


U!--55  =  A  ♦  BV|ne 


Data  Reduction.  The  two  assumptions  made  when  reducing  the  data 


are  (1)  the  boundary  layer  is  two-dimensional  and  (2)  the  instantaneous 


velocities  seen  by  the  two  parallel  wires  are  equal.  The  first  assumption 


is  needed  since  the  third  component  of  velocity  (w)  cannot  be  determined. 


If  the  probe  is  aligned  with  the  flow,  however,  and  if  the  boundary  layer  is 


two-dimensional,  the  w-component  makes  only  second  order  contributions  to 


sensor  response  and  may  be  safely  neglected.  If  the  second  assumption 


holds,  the  velocity  term  (Ueff)  may  be  eliminated  from  the  response 


equations  for  the  two  parallel  wires,  resulting  in  an  equation  solely  in 


terms  of  the  ambient  temperature,  T..  The  ambient  temperature  may  be 


found  by  iteratively  solving  for  T»  using  the  Newton-Raphson  method. 


The  instantaneous  velocities  U  and  V  may  then  be  found  from  the  two 


outer  orthogonal  wires  using  the  relation  given  by  Champagne  (1967). 


Knowing  the  instantaneous  values  of  U,  V,  and  T,  the  rms  fluctuation 


quantities  (u',v',  and  t')  and  the  cross  correlations  (u'v',  u't',  and  v't') 


may  be  determined.  The  equations  used  are  given  by 


Q.  -  /a.2  _  Asl  _  iLRxll 

q  1  *  Vi  N(N-1) 


N-1  NCN-1 


where  q^=u,  v,  or  t,  and  N=number  of  samples. 


' v't". 
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N  V  -  “  •  •  *  -  ■’  O  O."  «  '  H  * 


1 1 


s 


Calibration  Facility.  A  schematic  of  the  calibration  facility  is 
shown  on  Figure  13.  The  calibration  tunnel  was  designed  and  built  by  Wilson 
(1970).  Temperature  control  was  added  as  part  of  this  project.  A  brief 
description  follows. 

Compressed  air  at  a  line  pressure  of  120  psi  passes  through  a  series  cf 
cleaning  and  regulating  equipment  before  reaching  the  air  heater.  The 
heater  is  a  coil  taken  from  a  household  hair  dryer  inserted  into  a 
plexiglass  tube.  A  115  VAC  wall  outlet  regulated  by  a  voltage  regulator 
(Variac)  powers  the  coil.  Overheating  is  prevented  by  a  thermal  switch 
which  opens  the  circuit  if  the  temperature  becomes  excessive.  The  heated 
air  then  passes  into  an  insulated  10  gallon  mixing  tank.  Since  the 
temperature  of  the  flow  varies  with  the  flow  rate  through  the  coil,  and 
since  an  isothermal  flow  is  needed,  the  total  flow  rate  through  the  system 
must  be  kept  constant.  This  may  be  achieved  by  keeping  the  pressure  in  the 
tank  constant.  The  vent  and  control  valves  shown  on  the  schematic  perform 
this  function.  The  control  valve  regulates  the  flow  to  the  calibration 
tunnel,  while  the  vent  valve  is  adjusted  until  the  pressure,  as  measured  by 
a  manometer,  returns  to  a  preset  value.  A  relief  valve  prevents  excessive 
pressure  buildup  in  the  tank. 

The  calibration  tunnel  consists  of  a  stilling  chamber  with  turbulence 
reducing  screens  and  a  honeycomb  flow  straightener ,  and  an  elliptical 
contraction  followed  by  a  calibration  chamber  and  a  nozzle.  The  flow 
exiting  the  nozzle  was  found  to  have  a  turbulence  intensity  of  0.31.  The 
flow  could  be  continuously  varied  from  0  m/s  to  as  high  as  130  m/s  with  the 
appropriate  choice  of  exit  nozzle.  The  temperature  was  found  to  be  very 
uniform  (.01°C),  and  could  be  varied  from  room  temperature  (~25°C)  to  over 


H0°C. 


Other  Instrumentation. 

Pressure  Transducer:  The  transducer  used  is  of  the  variable 
reluctance  type  (Vaiidyne  DP^5),  and  is  designed  for  low-range,  differential 
measurements  of  extremely  low  pressure  differences.  The  diaphram  used  has  a 
pressure  range  of  8.9  cm  (3.5  in.  H2O)  with  an  accuracy  of  0.5$  full  scale. 
The  response  was  found  to  be  very  linear  throughout  the  entire  operating 
range.  Calibration  of  the  transducer  was  performed  against  a 
mi cro-manometer  (Dwyer  Microtector)  with  an  accuracy  of  ±0.06  mm  H2O.  The 
calibration  curve  was  obtained  by  a  linear  fit  to  the  data. 

Carrier  Demodulator:  The  anolog  output  of  the  pressure 
transducer  is  sent  to  a  demodulator  (Vaiidyne  CD-15)  which  provides  a  stable 
DC  output.  The  gain  and  offset  on  the  demodulator  enabled  tailoring  of  the 
signal  to  maximize  digitizing  resolution. 

A/D  Converter:  Three  A/D  converters  were  used — an  HP  3^3? A 
system  voltmeter  and  two  Norland  Prowler  digital  oscilloscopes.  The  former 
is  a  3^  digit  successive  approximation  digitizer  capable  of  sampling  up  to 
3600  samples  per  second  (ASCII  mode).  The  latter  is  a  12  bit  digitizer 
capable  of  simultaneously  sampling  two  channels  of  data  at  speeds  up  to  100 
kHz.  Two  buffers  store  up  to  t)096  data  points  each.  Two  Prowlers  may  also 
be  hooked  together  as  master  and  slave,  enabling  four  channels  of  data  to  be 
taken  simultaneously.  Simple  processing  of  the  data  can  be  done  within  the 
Norland  using  the  built-in  math,  calculus,  and  signal  processing  routines, 
or  the  data  may  be  sent  via  the  IEEE-ii88  interface  bus  to  the  computer. 

Computer;  A  Hewlett  Packard  Series  200  Model  16  personal 
computer  with  512  K  memory  was  used  for  data  reduction  and  as  a  controller. 
The  computer  is  linked  with  a  dual  disk  drive  (HP  9121)  which  utilizes  3*5 
in.,  256  K,  single-sided  micro  diskettes,  and  with  an  inkjet  printer  (HP 
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Thinkjet).  High  level  HP  Basic  is  the  programming  language  used.  A 
co-processor  (Infotek  FP210)  and  a  Basic  Compiler  (Infotek  BC203)  were  aaded 
to  speed  data  reduction.  A  5-10  times  increase  in  speed  was  realized.  An 
IEEE-488  interface  bus  allowed  communication  with  other  instruments. 

Anemometer:  A  four-channel  constant-temperature  anemometer 
(TSI  IFA-100)  was  used  to  drive  the  hot-wires.  The  anemometer  features  a 
built-in  mi croprocessor  with  non-volatile  memory  which  monitors  and  stores 
set-up  parameters  for  future  use.  The  high  frequency  response  and  high 
signal-to-noise  ratio  gives  accurate  measurements  even  in  low  turbulence 
intensity  flows.  A  built-in  square  wave  generator  enables  the  frequency 
response  to  be  determined.  Also,  the  built-in  signal  conditioners  are  used 
to  tailor  the  output  signal,  maximizing  the  digitizing  resolution  of  the  A/D 
converters . 

Calibration  of  the  Probe.  The  probe  must  be  calibrated  as  a 
function  of  both  velocity  and  temperature.  It  was  decided  to  calibrate  the 
probe  for  velocity  at  one  temperature,  then  change  temperatures ,  thus 
obtaining  a  calibration  curve  of  anemometer  voltage  vs.  velocity  for  each 
temperature  setting. 

A  schematic  of  the  set-up  used  is  shown  on  Figure  14.  The  probe  is 
driven  by  three  of  the  four  bridges  on  the  anemometer.  The  signal 
conditioners  offset  and  amplify  the  voltages,  maximizing  digitizing 
resolution.  Two  Norland  Prowlers,  wired  as  master  and  slave,  simultaneously 
sample  and  digitize  the  anemometer  signals.  Two  thousand  data  points  over 
eight  seconds  are  typically  taken.  The  math  processing  capability  of  the 
Prowlers  is  used  to  internally  average  the  signals,  after  which  the  averaged 
values  are  sent  to  the  computer  via  the  IEEE-488  data  link.  In  the 
meantime,  the  HP  34 37 A  A/D  converter  is  used  to  sample  the  pressure 


transducer  signal.  The  data  is  sent  to  the  computer  for  averaging.  Three 
thousand  points  taken  over  8  seconds  is  typical.  The  pressure  transducer 
signal  is  used  to  find  the  velocity  at  the  exit  of  the  calibration  jet 
through  a  calibration  equation  relating  the  transducer  signal  to  the 
pressure  drop  at  the  jet  exit.  The  jet  temperature  is  measured  with  a 
thermocouple.  The  Fluke  8840A  voltmeter  samples  the  thermocouple  voltage 
and  sends  the  data  to  the  computer  on  request.  Twenty  data  points  were 
taken  per  calibration  curve,  with  more  data  points  taken  at  lower  velocities 
where  the  calibration  curve  is  steep.  The  two  outer  wires  were  calibrated 
at  an  overheat  of  1.9,  while  the  center  wire  was  calibrated  at  an  overheat 
of  1.035. 

The  above  setup  enables  calibration  curves  to  be  found  quickly  and 
accurately.  The  process  is  highly  automated;  the  computer  controls  all  Cats 
acquisition,  with  the  experimenter  having  only  to  change  the  jet  velocity 
and  temperature. 

In  determining  the  calibration  coei f icients ,  the  response  equation 
given  by  equation  (2)  was  assumed  for  each  jet  temperature.  The  velocity 
Uc  is  not  the  jet  velocity,  but  the  effective  cooling  velocity  given  by 
Wilson  (1970) 

Uef f  *=  U[(1  +  kJ)0<ssin9]  (7) 

where  0  is  the  angle  between  the  wire  and  the  main  flow  direction.  The 
factor  k  is  introduced  to  account  for  cooling  due  to  the  velocity  tangential 
to  the  sensor  as  discussed  earlier.  The  wires  used  in  this  probe  had  1/d  = 
200.  The  coefficients  A  and  B  in  equation  (2)  are  obtained  using  a  least 
squares  fit  to  the  data  for  two  Jet  temperatures,  typically  25°C  and  30°C. 
The  variation  in  A  and  B  is  then  modelled  using  a  linear  fit  and  introduced 
into  equation  (2)  to  get  the  final  response  equation  given  by  equation  (3). 


Qualification  of  the  Probe.  Qualification  of  the  probe  took  place 
in  a  flat-plate  turbulent  boundary  layer  with  a  momentum  thickness  Reynolds 
number  of  Re0  -  2670  and  a  uniform  wall  heat  flux  boundary  condition 
(q"«226  W/m2).  The  probe  was  traversed  across  the  boundary  layer  and 
measurements  of  u'v',  t',  u't',  and  v't'  were  made.  These  quantities, 
normalized  on  the  free-stream  velocity  (Ux)  and  the  wall  to  free-stream 
temperature  difference  (Tw  -  TB),  were  compared  with  the  boundary  layer 
data  of  Blair  and  Bennett  ( 1 98-4 )  and  Gibson  and  Verriopoulos  (198*0.  All 
profiles  agreed  within  the  uncertainty  of  the  measurement  (estimated  to  be 
15%).  As  a  further  check  on  the  validity  of  the  data,  profiles  of  apparant 
t'  (the  t'  measured  by  the  probe  in  an  isothermal  boundary  layer)  were 
measured.  Apparant  t*  arises  due  to  the  two  parallel  wires  not  seeing  the 
same  velocity.  Apparant  t'  values  on  the  order  of  the  actual  t'  values 
would  indicate  a  problem  with  the  probe.  The  results  are  shown  on  Figure  15. 
It  may  be  seen  that  the  apparant  t*  values  are  5%  to  10t  of  the  actual  t' 
values  except  very  near  the  wall.  The  influence  on  the  actual  data  is  ever, 
smaller  since  the  apparent  t'  affects  the  actual  data  in  a  root-sum-square 
manner.  Interestingly  enough,  turbulence  intensity  is  seen  to  have  little 
effect  on  apparant  t’.  This  suggests  that  distance  from  the  wall  rather 
than  turbulence  intensity  is  the  controlling  parameter. 

Results  and  Discussion.  Measurements  were  made  in  the  flow 
discussed  in  section  A  of  this  report.  The  station  locations  are  shown  on 
F igure  1 . 

Low  TI  Case  (TI-0.fe5f):  Shear  stress  profiles  are  shown  on 
Figure  16.  The  profiles  are  seen  to  respond  quickly  to  curvature.  A 
dramatic  reduction  at  the  onset  of  curvature  is  observed,  after  which  the 
profiles  slowly  evolve  to  what  appears  to  be  an  asymptotic  state  by  station 
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6.  Curvature  is  seen  to  reduce  the  shear  stress  by  about  45?  below  the  flat 
wall  values  throughout  the  curve.  Recovery  of  shear  stress  is  slow,  with 
recovery  seemingly  complete  in  the  wake  by  station  10.  You,  et  al.  (1986) 
have  shown  that  the  near  wall  values  of  u'v'  recover  very  slowly.  All  the 
trends  are  consistent  with  the  observations  of  Gillis  and  Johnston  (1963) 
and  You,  et  al .  (1986). 

Profiles  of  t'  are  plotted  on  Figure  17.  The  flat  upstream  values 
(station  1)  are  seen  to  agree  with  the  measurements  of  Blair  and  Bennett 
(1984)  and  Gibson  and  Verriopoulos  (1984)  within  the  uncertainty  of  the 
measurement.  The  profiles  within  the  curve  assume  an  asymptotic  shape  for 
y/6  <  0.5.  No  evolution  of  the  profiles,  as  was  observed  for  the  shear 
stress,  is  evident.  Recovery  has  a  dramatic  effect  on  t'.  Profiles  of  t’ 
are  seen  to  overshoot  the  station  1  profiles.  It  is  believed  that  the 
profile  shape  eventually  returns  to  that  of  station  1 .  The  recovery  length 
was  too  short  to  observe  this  in  the  present  experiment,  however. 

The  effects  of  curvature  and  recovery  on  the  streanwise  turbulent  heat 
flux  u't'  are  shown  on  Figure  18.  The  station  1  values  agree  very  well  with 
the  data  of  Gibson  and  Verriopoulos  (1984).  The  effect  of  curvature  on  u't' 
is  dramatic.  The  profiles  in  the  curve  "snap"  into  an  asymptotic  shape  for 
y/6  <  0.5.  Slow  evolution  of  the  profiles  for  y/ <_  >  .05  may  be  observed. 
Recovery  is  seen  to  cause  an  overshoot  of  u't'  beyond  the  station  1  values, 
although  the  profile  is  expected  to  return  to  the  station  1  values  given  a 
long  enough  recovery  length. 

Profiles  of  the  apparent  turbulent  heat  flux  v't'  are  shown  on  Figure 
19.  The  effects  of  curvature  and  recovery  on  v't'  are  very  similar  to  those 
observed  for  u't' — a  dramatic  reduction  in  the  curve  to  an  asymptotic  shape, 
followed  by  an  overshoot  in  the  recovery.  Values  of  v't'  as  calculated  from 


the  wall  heat  flux  measurements  are  also  labelled  on  Figure  19  (at  y/6  = 
0.0).  The  wall  values  and  the  profile  measurements  agree  very  well  at 
station  1  and  in  the  curve.  Although  little  data  is  available  near  the  wall 
in  the  recovery  section,  the  values  of  v't'  seem  to  rise  above  the  wall  heat 
flux  values,  especially  for  station  10.  This  may  perhaps  be  seen  more 
clearly  on  Figure  20  where  v't'  has  been  normalized  on  wall  heat  flux  (q") 
rather  than  the  potential  velocity  at  the  wall  (UpW)  and  the  wall-to- 
freestream  temperature  difference  (Tw  -  T»).  The  profiles  at  station  1 
and  within  the  curve  approach  unity  near  the  wall,  while  the  profiles  in  the 
recovery  section  seem  to  overshoot.  It  is  believed  that  the  profiles 
overshoot,  then  turn  down  and  approach  unity  as  the  wall  is  approached. 

This  raises  the  question  of  whether  it  is  physically  realistic  for  the 
cross-stream  heat  flux,  v't',  to  increase  beyond  the  wall  heat  flux  when  no 
reversal  in  the  sign  of  v  and  t  is  observed.  A  possible  explanation  may  be 
found  in  the  time-averaged  energy  equation.  The  applicable  equation, 
neglecting  the  molecular  conduction  term,  o9t/3y,  is  given  by 


-at  -at  a(v't') 

u —  +  v —  =  -  - 

ax  ay  ay 


(8) 


If  v't'  is  to  increase  above  the  wall  heat  flux,  the  term  on  the  r.h.s.  of 
the  above  equation  must  be  negative.  The  term  v3t/ay  is  expected  to  be 
negative  since  v  >  0  (boundary  layer  is  growing),  and  3t/3y  <  0  (wall  is 
heated).  Since  u  is  always  positive,  the  condition  that  must  be  satisfied 
for  a  positive  gradient  in  v't'  is 


particular  y.  This  condition  could  occur  in  the  recovery  section.  You,  et 
al.  (1986)  and  Simon,  et  al.  (1980)  have  shown  that  the  Stanton  number 


increases  in  the  recovery,  implying  an  increase  in  the  heat  transfer 
coefficient,  h.  An  increase  in  h  with  x  implies  a  decrease  in  the 
temperature  near  the  wall,  i.e.,  3t/3x  <  0.  Physically,  eddy  diffusion  of 
heat  in  the  cross-stream  direction  is  increasing  at  the  expense  of 
convection  of  heat  in  the  streamwise  direction. 

Profiles  of  turbulent  Prandtl  number,  Prt,  are  shown  on  Figure  21.  The 
data  show  some  unavoidable  scatter,  but  it  is  seen  that  Pr^  decreases  in  the 
curve  below  the  station  1  values.  This  is  inconsistent  with  the 
observations  of  Gibson  and  Verriopoulos  (196^)  who  stated  that  Prt  rises 
within  the  curve  for  y/6  <  0.M  (though  this  was  not  conclusive  in  the  data 
due  to  scatter)  and  with  the  observations  of  You,  et  al .  (1986)  and  Simon, 
et  al .  (I960)  who  saw  Prt  rise  in  the  log-linear  region.  The  reason  for 
this  inconsistency  is  not  presently  known  although  it  should  be  noted  that 
the  second  and  third  references  present  Prt  values  deduced  from  the  mean 
temperature  and  velocity  profiles  which  apply  to  the  turbulent  core  whereas 
the  present  values  and  the  Gibson  and  Verriopoulos  values  were  measured  in 
the  wake. 

High  TI  Case  (TI-2.0E) : 

The  high  turbulence  case  was  obtained  by  inserting  a  grid  at  the 
entrance  to  the  contraction  section  as  discussed  in  section  A.  Merging  of 
the  boundary  layers  occurred  at  station  6  and  beyond  due  to  the  rapid  growth 
of  the  boundary  layer  on  the  unstable  concave  side  for  the  high  TI  case. 

The  data  taken  at  station  6  and  beyond  were  therefore  disqualified. 

Shear  strees  profiles  are  presented  on  Figure  22.  The  higher 
turbulence  intensity  has  the  effect  of  increasing  the  shear  stress 
throughout  the  boundary  layer.  This  is  expected.  The  profiles  immediately 
assume  an  asymptotic  shape  within  the  curve  for  y/6  <  0.5,  unlike  the  low  TI 


case  where  an  evolution  of  profiles  was  observed.  The  asymptotic  profile 
for  the  low  TI  case  is  also  plotted  on  Figure  22,  and  is  seen  to  agree  well 
with  the  high  TI  data  in  the  curve.  A  reversal  in  the  sign  of  the  shear 
stress  is  seen  to  occur  in  the  outer  part  of  the  boundary  layer,  along  with 
a  slow  relaxation.  This  reversal  in  sign  has  also  been  observed  by  Gill  is 
and  Johnston  (1983)  for  more  strongly  curved  flows.  The  cause  of  this 
reversal  is  not  yet  understood,  bu*-  is  thought  to  be  linked  to  the  sudden 
drop  in  the  production  of  shear  stress  in  the  curve.  This  is  presently 
under  study. 

The  effect  of  turbulence  intensity  on  t*  is  shown  on  Figure  23.  The 
flat  plate  values  are  seen  to  be  greater  than  the  low  turbulence  values  by 
approximately  60$.  This  increase  is  completely  eliminated  by  station  4, 
however,  where  tne  values  agree  with  the  asymptotic  profile  of  the  low  TI 
case.  This  dramatically  illustrates  the  dominance  of  curvature  over 
turbulence  intensity. 

The  dominance  of  curvature  is  again  illustrated  by  the  u't'  profiles 
shown  on  Figure  24,  Although  the  flat  plate  values  for  the  high  TI  case  are 
150$  higher  than  the  corresponding  values  for  the  low  TI  case,  a  profile 
agreeing  with  the  asymptotic  profile  for  the  low  TI  case  is  achieved  by 
station  4.  The  effect  of  curvature  is  remarkable. 

Profiles  of  v't'  are  given  on  Figure  25.  The  usual  trends  are  observed. 
Values  of  v't’  higher  than  the  corresponding  low  TI  values  are  seen  at 
station  1,  while  an  asymptotic  profile  agreeing  with  the  asymptotic  profile 
for  the  low  TI  case  is  achieved  in  the  curve,  v't'  is  seen  to  react  faster 
than  u't'  and  t',  however,  achieving  an  asymptotic  shape  by  station  3 
instead  of  station  4.  The  profiles  are  seen  to  agree  well  with  the  wall 
heat  flux  measurements.  Profiles  of  v’t'  normalized  on  the  wall  heat  flux 
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are  shown  on  Figure  26.  The  station  1  values  are  expected  to  level  out  and 
approach  unity  near  the  wall. 

Turbulent  Prandtl  number  profiles  are  shown  on  Figure  27.  A  slight 
decrease  in  Prt  within  the  curve  near  the  wall  is  seen.  It  is  interesting 
that  in  this  profile  the  values  of  Prt  become  negative  in  the  wake  due  to 
the  sign  reversal  in  shear  stress,  as  discussed  earlier.  This  strong 
variation  in  Prt  indicates  that  Reynolds  analogy  is  not  valid. 

Conclusions .  The  effects  of  curvature,  recovery  and  turbulence 
intensity  on  turbulenc  boundary  layer  heat  transfer  were  studied.  The  main 
conclusions  are: 

(1)  Profiles  of  u'v',  t',  u't',  and  v't'  are  reduced  in  the  curve 
with  an  asymptotic  profile  being  achieved  by  station  i4. 
Recovery  in  the  wake  region  occurs  quickly,  with  the  profiles 
often  overshooting  the  flat  upstream  (station  1)  values. 

(2)  A  reversal  of  sign  in  the  shear  stress  is  seen  for  the  high  TI 
case  (TI=2.0?)  for  y/6  >  0.5.  This  is  consistent  with  the 
observations  of  other  researchers. 

(3)  Curvature  is  seen  to  decrease  the  turbulent  Prandtl  number  in 
the  wake  region.  Negative  turDulent  Prandtl  numbers  were  seen 
in  the  wake  for  the  high  TI  case  and  Reynolds  analogy  is 
disqualified. 

(*0  Curvature  effects  dominate  turbulence  intensity  effects  for 
the  cases  studied. 

3.  Expected  Activities 

This  probe  development  has  been  completed.  More  fine  tuning  will,  of 
course,  be  done  as  the  probe  is  used  in  subsequent  experiments. 

D.  Transport  in  a  Turbulent  Boundary  Layer  with  Concave  Curvature 

1 .  Objectives 

This  study  investigates  the  influence  of  streamwise  concave  curvature 
on  momentum  and  heat  transport  in  a  fully  turbulent  boundary  layer. 


Parameters  such  as  freestream  turbulence  intensity  and  accelerat ion ,  which 
affect  the  flow  on  gas  turbine  blades,  will  be  included. 

This  study  is  directed  toward  a  better  understanding  of  the  flow  and 
heat  transfer  physics  with  emphasis  on  the  recovery  of  the  flow  following 
concave  curvature.  In  this  region,  flow  visualization  will  be  carried  out. 
Also,  the  existing  database  of  turbulence  and  heat  transfer  quantities  will 
be  extended  to  investigate  the  effect  of  different  radii  of  concave 
curvature  including  the  recovery  behavior  of  momentum  and  heat  transport. 

2.  Progress  and  Accomplishments 

A  fan  speed  unit  for  the  present  wind  tunnel  is  currently  being 
installed.  This  will  allow  calibrating  hot  wire  probes  for  velocity  and 
temperature  measurements  in-si tu  by  changing  the  speed  of  the  wind  tunnel. 
Also,  experimental  data  can  be  taken  for  diffei  ent  free-stream  velocities. 
The  speed  of  the  wind  tunnel  will  be  directly  controlled  from  a  computer. 

Literature  on  previous  studies  including  basic  studies  and  reports  on 
experimental  techniques  are  being  reviewed. 

3.  Expected  Activities 

An  experimental  setup  which  allows  incorporating  the  effect  of  flow 
acceleration  in  curved  ducts  of  various  radii  of  curvature  is  being  designed 
which  employs  bendable  duct  walls.  Merging  of  the  boundary  layers  at  the 
inner  (convex)  and  the  outer  (concave)  wall  will  be  prevented  by  applying 
suction  at  strategic  locations. 

Barlow  (1985)  carried  out  a  visualization  study  of  the  flow  over  a 
concave  surface  in  a  water  tunnel.  Roll  cells  which  develop  due  to  the 
inertial  instability  mechanism  were  visualized  in  the  curve.  The  present 
study  will  continue  this  work  by  investigating  the  behavior  of  the  roll 
cells  on  the  flat  recovery  section  downstream  of  the  curve.  Appropriate 
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flow  visualization  techniques  will  be  chosen.  A  liquid  crystal  technique 
as  used  by  Simmonich  and  Moffat  (1982)  allows  visualization  of  the  roll 
cells  if  they  are  fixed  in  position  by  vortex  generators. 

A  liquid  crystal  processing  method,  as  used  by  Hi ppensteele ,  Russell 
and  Torres  (1985)  for  turbine  blades,  will  be  considered  for  measuring  the 
local  heat  transfer  coefficient.  This  method  uses  the  isotherms 
corresponding  to  the  yellow  color  of  a  liquid  crystal  sheet  on  a  uniform 
heat  flux  wall  to  determine  the  local  heat  transfer  coefficient. 

Also,  a  technique  using  naphthalene  dissolved  in  methyl  alcohol  and 
sprayed  on  the  wall  (McCormack,  Welke  and  Kelleher — 1970)  may  be  applied. 
Smoke  techniques  are  being  considered  to  observe  the  moving  roll  cells  in 
natural  flow  (without  vortex  generators). 

E.  Film  Cooling  on  Curved  Surfaces 

This  study  is  concerned  with  the  film  cooling  on  concave  and  convex 
surfaces  of  constant  radius  of  curvature.  It  is  desired  to  observe  the 
effect  of  the  ratio  of  radius  of  curvature  to  half  the  injection  hole 
diameter  (2r/D  —  strength  of  curvature)  on  film  cooling  effectiveness.  A 
wind  tunnel  and  test  section  (discussed  in  previous  progress  reports)  has 
been  constructed  for  this  purpose.  The  test  section  has  concave  and  convex 
walls  with  three  sets  of  replaceable  injection  holes.  Each  set  consists  of 
a  row  of  holes  inclined  at  35  degrees  to  the  mainflow.  Hole  diameters  were 
such  that  2r/D  on  the  concave  surface  was  either  -57,  -89,  or  -111.  On  the 
convex  surface  they  were  61,  9^,  or  115. 

Effectiveness  measurements  have  been  taken  for  all  sets  of  injection 
geom.ter ies  on  both  surfaces  for  density  ratios  of  2.0  and  0.95.  Figure  27 
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shows  the  lateral-averaged  effectiveness  far  downstream  of  injection 
( x / D= -4 0 )  with  a  density  ratio  of  0.95. 

Stronger  convex  curvature  is  found  to  enhance  effectiveness  at  low 
blowing  rates,  M,  while  at  high  blowing  rates,  strength  of  curvature  seems 
to  diminish  in  importance.  A  blowing  rate  is  found  where  an  increase  in 
itsmagnitude  causes  a  decrease  in  ef fecti veness .  The  centrifugal  force  of 
the  jet  at  high  M  causes  the  jet  to  lift  off  the  surface.  Further  increase 
in  blowing  rate  causes  the  jets  to  merge  and  blockage  of  the  mainflow  from 
the  wall  results.  An  increase  in  effectiveness  is,  therefore,  encountered. 

On  the  concave  wall,  at  low  blowing  rates,  effectiveness  is  virtually 
the  same  for  all  strengths  of  curvature,  while  at  higher  blowing  rates, 
increased  curvature  strength  raises  effectiveness.  An  increase  in  blowing 
rate,  on  the  concave  surface  always  causes  an  increase  in  effectiveness;  the 
centrifugal  force  of  the  jets  impinges  on  the  concave  wall. 

The  flat  plate  study  of  Pedersen  (1972)  is  included  in  Figure  26  for 
comparison . 

F.  End-Wall  Heat  Transfer 

1 .  Introduction 

This  work,  employing  the  heat-mass  transfer  analogy  using  naphthalene 
sublimation,  is  an  investigation  of  mass  (heat  transfer)  on  the  end-wall. 

The  mass  transfer  technique  offers  much  greater  detail  of  the  local  transfer 
coefficient  than  if  heat  transfer  measurements  were  taken  directly. 

2.  Progress  and  Accomplishments 

In  the  previous  progress  report,  results  were  given  for  the  case  of  Re 
«  1 . M  x  1 05  (based  on  chord  length  of  the  cascade)  and  boundary  layer 


thickness  of  6  *  0.592  in.  This  is  referred  to  as  case  1.  In  this  update, 
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contour  plots  are  given  for  two  additional  test  cases  plus  some  revealing 
one-dimensional  heat  transfer  plots  for  the  endwall  region. 

Test  cases  2  and  3,  shown  in  Figures  29  and  32,  respectively,  were 
conducted  under  the  fallowing  conditions:  case  2  -  Re  =  1.42  x  10^,  6  = 
1.147  in.  and  case  3  *  Re  =  8.86  x  1014,  6  =  0.461  in.  A  comparison  between 
cases  1  and  2  shows  the  effect  of  an  increased  boundary  layer  thickness 
while  comparing  cases  1  and  3  shows  the  Re  dependence.  Note  that  the  data 
is  presented  as  a  ratio  of  the  local  mass  transfer  on  the  endwall  to  the 
mass  transfer  that  would  be  found  if  no  turbine  blades  were  present  (pure 
flat  plate  situation).  Although  the  contour  maps  indicate  differences  in 
heat  transfer  between  the  three  cases  of  up  to  ~  25?  for  certain  regions, 
the  data  averaged  across  the  passage  shows  virtually  no  differences  (see 
Figure  31 ) . 

Two  one-dimensional  plots  are  also  informative.  First,  Figure  32  shows 
the  variation  in  heat  transfer  at  locations  upstream  from  the  leading  edge 
of  a  turbine  blade.  Note  the  two  local  maxima  in  the  graph  showing  the 
effects  of  the  inner  vortex  (higher  maximum)  and  the  main  horseshoe  vortex 
(smaller  maximum).  Figure  33  illustrates  how  the  heat  transfer  coefficient 
changes  across  the  passage.  The  extreme  drop  in  heat  transfer  occurs  where 
the  endwall  flow  separates  to  climb  the  suction  surface  and  become  entrained 
in  the  passage  vortex,  while  the  maximum  near  the  suction  surface  is  due  to 
the  formation  of  a  suction  side  corner  vortex. 

3.  Expected  Activities 

The  writing  of  a  formal  paper  on  the  endwall  phase  of  this  project  is 
currently  underway  while  the  development  of  equipment  to  measure  mass 
transfer  on  turbine  blades  continues. 


G.  Near-End-Wall  Heat  Transfer 
1 .  Objectives 

The  present  study  is  to  investigate  the  performance  of  film  cooling  on 
the  convex  surface  of  the  turbine  blade  through  two  rows  of  injection  holes 
in  the  neai — end-wall  region.  Several  parameters  which  do  not  occur  at  the 
two-dimensional  flow  region  do  affect  the  performance  of  film  cooling. 

These  parameters  include:  lower  free  stream  velocities  and  higher  turbulent 
intensity  inside  the  boundary  layer,  the  passage  vortex  traveling  from  the 
concave  side  of  the  adjacent  blade,  horseshoe  vortex  generated  from  the 
leading  edge,  deflected  mainstream  direction,  and  the  plugged  injection  hole 
of  the  second  row.  The  passage  vortex,  which  plays  a  dominant  role,  sweeps 
the  coolant  from  the  convex  surface  and  leaves  a  triangular  region  without 
film  protection. 

Test  facilities  and  detailed  arrangements  of  injection  holes  and 
sampling  taps  are  shown  in  Figures  3^  and  35,  respectively.  The  injection 
hole,  half  blocked  by  the  end  wall,  is  plugged  to  assure  the  uniformity  of 
secondary  jet  through  each  injection  hole.  Results  with  lower  density 
ratio,  R=0.96,  were  reported  last  year.  Results  on  the  convex  surface  with 
higher  density  ratio,  R-2.0,  are  presented  herein. 

2.  Progress  and  Accomplishments 

Presently,  the  approaching  velocity  which  is  measured  at  12.9  cm  from 
the  leading  edge  of  the  turbine  blade  is  9.8  m/s.  The  boundary  layer 
thickness,  699,  at  the  same  location  is  2  cm  and  the  turbulence  intensity, 
Tu,  of  the  freestream  is  1.2$.  The  Reynolds  number,  Re*,  is  based  on  the 
freestream  velocity,  U„>=17.8  m/s,  at  the  injection  hole  of  the  first  row 
and  the  diameter  of  the  injection  hole.  Figures  36  through  38  show  the 
values  of  the  average  film  cooling  effectiveness  at  blowing  rates  of  M=0.5, 


1.02,  and  1.^7,  respectively.  The  solid  symbols  on  the  right  axis  are 
values  of  average  effectiveness  measured  in  the  two-dimensional  flow  region 
at  midspan  of  the  blaae.  Tne  open  and  solid  circular  symbols  are  mapped 
into  the  convex  surface  of  the  turbine  blade  and  form  the  boundaries  of 
region  B  and  region  C  in  Figures  39  and  41.  The  distributions  of  average 
effectiveness  at  density  ratios  of  R=2.0  and  0.96  are  shown  in  Figures  39 
and  HO,  respectively. 

Region  C,  where  coolant  is  swept  from  the  surface  by  the  effect  of  the 
passage  vortex,  is  the  unprotected  region.  In  region  B,  the  values  of 
average  effectiveness  at  each  downstream  location  decrease  linearly  to  zero 
as  the  elevation,  H,  moves  from  the  upper  boundary  toward  the  lower  boundary. 
As  expected,  the  area  of  the  region  C  slightly  decreases  with  increasing 
blowing  rates  and  the  upper  boundary  of  the  region  B  stays  the  same  as  the 
blowing  rate  changes.  Comparison  of  the  distribution  of  average 
effectiveness  between  two  different  density  ratios  R=2.0  and  R=0.96  is  shown 
in  Figure  Hi.  Boundaries  of  region  B  move  slightly  toward  the  end  wall. 
Comparison  of  values  of  the  average  effectiveness  is  shown  in  Figure  H2. 

3.  Expected  Activities  Next  Contract  Year 

Complete  results  on  the  concave  surface  with  higher  density  ratios  will 
be  presented  next  year. 
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Table  1.  Descriptors  of  CASES  1,  2  &  3 
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Figure  1.  Plan  View  of  the  Curved  Boundary  Laver  Facility 
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Figure  12.  Schematic  of  probe 


Figure  13.  Schematic  of  calibration  facility 
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Figure  17.  t*  profiles  (TI  =  0.65%) 


Figure  18.  u't’  profiles  (TI  =  0.65%) 
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Figure  31.  Comparison  of  average  data.  Cases 
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Figure  32.  St  /St  vs.  distance  from  leading  edge 
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Figure  35.  Detail  of  test  blade  showing  measurement  regions 
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Figure  37.  Average  film  cooling  effectiveness  on  convex  surface  at 


Distribution  ranges  of  average  film  cooling  effectiveness 
on  the  convex  surface  at  near-end-wall  region 
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Figure  42.  Comparison  of  values  of  average  film  cooling  effectiveness 
between  two  different  density  ratios 
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